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ABSTRACT 

We present the results from multi-frequency VLBA observations from 327 MHz to 8.4 
GHz of the gigahertz-peaked spectrum radio source PKS 1518+047 (4C 04.51) aimed 
at studying the spectral index distribution across the source. Further multi-frequency 
archival VLA data were analysed to constrain the spectral shape of the whole source. 
The pc-scale resolution provided by the VLBA data allows us to resolve the source 
structure in several sub-components. The analysis of their synchrotron spectra showed 
that the source components have steep spectral indices, suggesting that no supply/re- 
acceleration of fresh particles is currently taking place in any region of the source. 
By assuming the equipartition magnetic field of 4 mG, we found that only electrons 
with 7 600, are still contributing to the radio spectrum, while electrons with higher 
energies have been almost completed depleted. The source radiative lifetime we derived 
is 2700±600 years. Considering the best fit to the overall spectrum, we find that the 
time in which the nucleus has not been active represents almost 20% of the whole 
source lifetime, indicating that the source was 2150±500 years old when the radio 
emission switched off. 
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1 INTRODUCTION 

The radio emission of extragalactic powerful radio sources 
is due to synchrotron radiation from relativistic particles 
with a power-law energy distribution. They are produced 
in the "central engine", namely the active galactic nucleus 
(AGN), and reaccelerated in the hot spot, that is the region 
where the particles, channelled through the jets, interact 
with the external medium. The energy distribution of the 
plasma is described by a power-law: N(E) — N E~ S that 
results in a power-law radio spectrum S v cx v~ a , with a 
moderate steepening S v oc v ~ a + - 5 caused by energy losses 
in the form of radio emission. The spectral steepening 
occurs at the break frequency, Vh, whi ch is related to th e 
age of relativistic electrons (see e.g. iPacholczvk I Il970h . 
This model, known as the continuous injection (CI) model, 
requires a continuous injection of power-law distributed 
electrons into a volume permeated by a constant magnetic 
field. The model has been applied to the in terpretation of 
the total spectra of the radio sources (see e.g. lMurgia et alTI 
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1999), assuming that the emission of the lobes (supposed 
to have a constant and uniform magnetic field) dominates 
over that of core, jets and/or hot spots. 
It is nowadays clear that powerful radio sources represent a 
small fraction of the population of the active galactic nuclei 
associated with elliptical galaxies, implying that the radio 
activity is a transient phase in the life of these systems. The 
onset of radio activity is currently thought to be related to 
merger or accretion events occurring in the host galaxy. 
The evolutionary stage of each powerful radio source is 
indicated by its linear size. Intrinsically compact radio 
sources with linear size LS ^ 1 kpc, are therefore inter- 
preted as young radio objects in which the radio emission 
originated a few tho usand years ago (e.g. iFanti et al. I ll995l : 
ISnellen et aL~ll200Ch . These sources are characterised by a 
rising synchrotron spectrum peaking at frequencies around 
a few gigahertz and are known as gigahertz-peaked spec- 
trum (GPS) objects. When imaged with high resolution, 
they often display a two-sided morphology which looks like 
a scaled-down version of the cla ssical and large (u p to a 



few Mpc) double radio sources dPhillips fc Mutel I fl982h 
For these sources, kinematic (|Polatidis fc Conway 1 12003) 
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and radiative i|Murgiall2003h studies provided ages of the 
order of 10 3 and 10 4 years, supporting the idea that they 
are young radio sources whose fate is probably to become 
classical doubles. However, it has been claimed (jAlexander I 
2000:; iMarecki et al. 1 120061 ) that a fraction of young radio 
sources would die in an early stage, never becoming the 
classical extended radio galaxies with linear sizes of a few 
hundred kpc and ages of the order of 10 7 - 10 8 yr. Support 
for this sc enario comes from a s tatistical study of GPS 
sources by iGugliucci et~al~l (|2005l ): they showed that the 
age distribution of the small radio sources they considered 
has a peak around 500 years. 

So far, only a few objects have been recognised as dy- 
ing radio sources. They are difficult to find owing to 
their extremely steep radio spectrum which makes them 
under-represented in flux-limited catalogues . Th e ob- 
jects 0809 +404 dKunert-Bairaszewska et al. I 120061 ) and 
1542+323 l|Kunert-Bairaszewska et al. I l2005h are possible 
examples of young radio sources that are fading. These 
sources, with an estimated age of 10 4 - 10 5 years, have been 
suggested as "faders" due to the lack of active structures, 
like cores and hot-spots, although their optically-thin radio 
spectra do not display the typical form expected for a fader. 

In this paper, we present results from multi- 
frequency VLBA and VLA data of the GPS radio 
source PKS 1518+047 (RA=15 h 21 m 14 s and DEC= 
04° 30' 2l" J2000), identified w ith a quasar at redshift 
z = 1.296 (jStickel fe Ruhr 1 1 1996] ). It is a powerful radio 
source (Log Pi.4GHz (W/Hz) = 28.53), with a linear size of 
1.28 kpc. This source wa s selected from the GPS sample of 
IStanghellini et al~1 |T998) on the basis of its steep spectrum 
Q i!4GHz = 1-2, uncommon in young radio sources. The goal 
of this paper is to determine by means of the analysis of the 
radio spectrum and its slope in the optically thin region, 
whether this source is in a phase where its radio activity 
has, perhaps temporarily, switched off. 

Throughout this paper, we assume the following 
cosmology: Ho = 71kms~ 1 Mpc~ , Qm = 0.27 and 
Qa ~ 0.73, in a flat Universe. At the redshift of the target 
1 = 8.436 kpc. The spectral index is defined as S(v) oc v~ a . 



2 RADIO DATA 

The target source was observed on March 5, 2008 (project 
code BD129) with the VLBA at 0.327, 0.611 (P band), 1.4 
and 1.6 GHz (L band) with a recording band width of 16 
MHz at 256 Mbps for 20 min in L band and 1 hr in P band. 
The correlation was performed at the VLBA correlator in 
Socorro. These observations have been complemented with 
archival VLBA data at 4.98 (C band) and 8.4 GHz (X 
band) carried out on March 28, 2001 (project code BS085). 
The data reduction was carried out with the NRAO AIPS 
package. 

During the observations, the gain corrections were found to 
have variations within 5% in L, C and X bands, and 10% 
in P band respectively. In P band, the antennas located 
in Kitt Peak, Mauna Kea and North Liberty had erratic 
system temperatures, probably due to local RFI, and their 
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Table 1. Multi-frequency flux density of PKS 1518+047 and of 
its northern and southern componen ts, a = archival VLA data ; b 
= VLBA data; c = WSRT data fromlStanghellini et al~l dl998h ; d 
= R ATAN data fr o mlStTanghellini et al. I jl998h : e = WSRT data 
from lXiang et al. I ||2006| ). 



data had to be flagged out completely. The resulting loss 
of both resolution and sensitivity did not affect the source 
structure and flux density. 

To constrain the spectral shape of the whole source, we 
analysed archival VLA data at 0.317, 1.365, 1.665, 4.535, 
4.985, 8.085, 8.465, 14.940 and 22.460 GHz, obtained on 
August 22, 1998 (project code AS637). The data reduction 
was carried out with the standard procedure implemented 
in the NRAO AIPS package. Uncertainties in the deter- 
mination of the flux density are dominated by amplitude 
calibration errors, which are found to be around 3% at all 
frequencies. 

VLA and VLBA final images were obtained after a number 
of phase-only self-calibration iterations. At 4.9 and 8.4 GHz, 
besides the full resolution VLBA image, we also produced 
a lower resolution image using the same Mi>-range, image 
sampling and restoring beam of the 1.6 GHz, in order to 
produce spectral index maps of the source. Images at the 
various frequencies were aligned using the task LGEOM 
by comparing the position of the compact and bright 
components at each frequency. 

Flux density and deconvolved angular sizes were measured 
by means of the task JMFIT which performs a Gaussian 
fit to the source components on the image plane, or, in the 
case of extended components, by TVSTAT which performs 
an aperture integration on a selected region on the image 
plane. 

The parameters derived in this way are reported in Table 
[1] About 8% and 20% of the total flux density is missing 
in our VLBA images at 4.9 and 8.4 GHz respectively if we 
consider the VLA measurements. It is likely that such a 
missing flux density is from a steep-spectrum and diffuse 
emission which is not sampled by the VLBA data due to 
the lack of appropriate short spacings. 
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3 RESULTS 



3.2 The spectral shape 



3.1 Source structure 

Multi-frequency VLBA observations with parsec-scale res- 
olution allow us to resolve the structure of PKS 1518+047 
into two main source components (Fig. [TJ, separated by 
135 mas (1.1 kpc), i n agr ee ment with pr e vious images 
by Dallacasa et al I (|l998l ). IXiang et all (|2002h . and 
IXiang et al. I ((2006). Both the northern and southern 
components have angular size of 23x15 mas 2 . 
At the highest frequencies, where we achieve the best 
resolution, the northern component is resolved into two 
sub-structures (labelled Nl and N2) separated by 11 mas 
(90 pc) and position angle of 11°. Nl accounts for 295 and 
119 mjy at 4.9 and 8.4 GHz, respectively, with angular size 
of 7.7x6.8 mas 2 . Component N2 is fainter than Nl and it 
accounts for 97 and 29 mjy at 4.9 and 8.4 GHz, respectively, 
with angular size of 7.4 x5 mas 2 . The southern complex 
is resolved into 4 compact regions located at 3, 7 and 16 
mas (25, 60 and 130 pc for S2, S3 and S4 respectively) 
with respect to SI, all with almost the same position angle 
of 40°. The resolution, not adequate to properly fit the 
several components of the southern lobe, does not allow us 
to reliably determine the observational parameters of each 
single sub-component. 

The spectral index in both the northern and southern 
complexes, computed considering the integrated component 
flux density, is steep, with a\'% = 1.5 ± 0.1 and = 1.0 ± 0.1, 
respectively. Errors on the spectral indices were calculated 
following the error propagation theory. The analysis of 
the spectral index distribution in the southern complex 
(Fig. [2]) shows a steepening of the spectral index going 
from SI inwards to S4, suggesting that SI is likely the last 
place where relativistic electrons coming from the source 
core were re-accelerated. In the northern lobe, the spectral 
index between 1.6 and 4.9 GHz is almost constant across 
the component, showing a small gradient going from Nl 
towards N2, while in the spectral index maps at higher 
frequencies it steepens from Nl to N2 (Fig. 0. The steep 
spectral indices showed by the northern and southern 
components suggest that no current particle acceleration 
across the source is taking place, indicating that active 
regions, like conventional jet-knots and hot spots, are no 
longer present. 

At 327 and 611 MHz, i.e. in the optically-thick part of 
the spectra, the low spatial resolution does not allow us 
to separate the individual sub-components. However, their 
spectral indices, obtained considering the integrated flux 
density, are 0327 ~ -0.6±0.2 and -0.9±0.2 for the northern 
and southern component respectively, that is much flatter 
than the value expected in the presence of "classical" syn- 
chrotron self-absorption from a homogeneous component, 
indicating that what we see is the superposition of the 
spectra of the various sub-components, each characterised 
by its own spectral peak occurring in a range of frequencies 
and then causing the broadening of the overall spectrum 
(i.e. the northern and southern components are far from be- 
ing homogeneously filled by magnetised relativistic plasma) . 



To understand the physical processes taking place in this 
source we fit the optically-thin part of the overall spectrum, 
as well as the spectra of the northern and southern com- 
ponents, assuming two different models. The first model 
assumes that fresh relativistic particles are continuously 
injected in the source (CI model), while in the second model 
the continuous supply of particles is over and the radio 
source is already in the relic phase (CI OFF model). 



The CI synchrotron model is described by three param- 



eters: 



the injection spectral index; 



ii) Vb, the break frequency; 

iii) norm, the flux normalization. 

As the injection of fresh particles stops (CI OFF model) , 
a second break, v\,2, appears at high frequencies, and beyond 
that the spectrum cuts off exponentially. This second break 
is related to the first according to: 

t \ 2 

(1) 

total source's age and toFF is 
since the i njection switched off 
(see e.g. iKomissarov fc Gubanov 1 1 19941 : ISlee et al. I l200ll : 



Vh2 = v h 

where t B is the 
the time elapsed 



I Parma et al. 2007). Indeed, compared to the basic CI 



model, the CI OFF model is characterized by one more free 
parameter: 

i) oiinj, the injection spectral index; 

ii) Vb, the lowest break frequency; 

iii) norm, the flux normalization; 

iv) ioFF/ts, i.e. the relic to total source age ratio. 

The spectral shapes of both models cannot be described 
by analytic equations and must be computed numerically 
(see e.g. Murgia et al. 1999 and Slee et al. 2001 for further 
details). 

In the study of the overall spectrum, in addition to 
the archival VLA data analysed in this paper, we con- 
sider also RATAN observ ations at 0.96, 3.9, and 11.1 GHz 
|Stanghellini et al. I B-9981 ), and WSRT observations at 5 
GHz ( Xiang et al. I [2006 s ) . to have a better frequency sam- 
pling. Due to the presence of another radio source located 
at about 1 arcminute from the target (Fig. [3]), we do not 
consider observations with resolution worse than ~1 . 
During the fitting procedure, a particular care was taken 
in choosing the most accurate injection spectral index. The 
peculiar shape of the radio spectrum of this source does not 
allow us to directly derive the injection spectral index from 
the optically-thin emission below the break, because it falls 
below the peak frequency where the spectrum is absorbed. 
For this reason, we fit the spectrum assuming various in- 
jection spectral indices, choosing the one that provides the 
best reduced chi-square (Fig. [4p). 

The best fit we obtain with the CI model implies a very 
steep injection spectral index a ln j = 1.1, that is reflected 
in an uncommon electron energy distribution 8 = 3.2, (Fig. 
13^). On the other hand, the CI OFF model provides a more 
common injection spectral index a- m j — 0.7 (Fig. |4p). Fur- 
thermore, the comparison between the reduced chi-square of 
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Figure 1. VLBA images at 0.327, 0.611, 1.6, 4.9 and 8.4 GHz of PKS 1518+047. On the images we show the observing frequency, the 
peak flux density and the first contour intensity (f.c.) which is 3 times the lc noise level measured on the image plane. Contours increase 
of a factor 2. The beam is plotted on the bottom left corner. 



the different models (Fig. [4}j) shows that the CI OFF model 
is more accurate in fitting the data. In the context of the CI 
OFF model, we derive a break frequency fbr = 2.4 GHz, and 
toFF/ts = 0.2. A similar result is found for the spectrum of 
the southern component. As in the overall spectrum, both 
CI and CI OFF models well reproduce the spectral shape 
(Fig. [5} , but the former implies an uncommonly steep injec- 
tion spectral index. 

A different result is found in the analysis of the northern 
lobe. In this case, we find that the CI OFF model with 
ai n j = 0.7 well reproduces the spectral shape, providing a 
break frequency — 0.8 GHz and ioFF/ts ~ 0.27 (Fig. 
[6p), while the CI model provides a worse chi-square even 
considering a steep injection spectral index (Fig. [6^.). 

The optically-thick part of the spectra is well modelled 
by an absorbed spectrum with athlck = —1.2 ± 0.1, that 
is different from the canonical -2.5 expected in the pres- 
ence of synchrotron self-absorption from a homogeneous 
component. This result, together with the source structure 
resolved in several sub-components, indicates that in the 
optically-thick regime the observed spectra are the super- 
position of the spectra of many components that cannot be 



resolved due to resolution limitation. 

The best fit parameters obtained with the CI and CI OFF 
models, and assuming synchrotron self-absorption (SSA) 
are reported in Table [2] 



3.3 Physical parameters 

Physical parameters for the radio source components were 
computed assu ming equipartition condition and using stan- 
dard formulae (|Pacholczvklll970l) . We considered particles 
with energy between 7 m i n =100 and 7 ma x =600. The high 
energy cut-off corresponds to the value for which the break 
frequency in the observer frame occurs at 2.4 GHz, as ob- 
tained from the best fit to the model (Section 3.2). Pro- 
ton and electron energy densities are assumed to be equal, 
and the spectral index is a = 0.7 (see Section 3.2). We as- 
sume that the volume V of the emitting regions is a prolate 
spheroid: 

V = £ ab 2 <t> 
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Figure 2. VLBA spectral index images of PKS 1518+047 between 1.6 and 4.9 GHz (panel a), between 4.9 and 8.4 GHz (panel 6), and 
between 1.6 and 8.4 GHz (panel c). On the images we show the observing frequency of the contours, the peak flux density and the first 
contour intensity (fx.) which is 3 times the la noise level measured on the image plane. Contours increase of a factor 2. The beam is 
plotted on the bottom left corner. The grey scale is shown by the wedge at the top of each image. 



where a and b are the spheroid major and minor axes, 
and <f) is the filling factor. In the case of the entire radio 
source, we consider a =140 mas, b =25 mas, and <j> — 1. 
With these assumptions, we obtain an equipartition mag- 
netic field H eq = 4 mG, and a minimum energy density 
itmin = 1.4 x 10~ 6 erg/cm 3 . In the case of component N and 
S we assume a = 23 mas, b =15 mas and cj>=l, and we infer 
a magnetic field of 7 mG, and a minimum energy density 
of =5xl0 -6 erg/cm 3 , in agreement with the values derived 
for the entire source structure. 



4 DISCUSSION 

In general the shape of the synchrotron spectrum of active 
radio sources is the result of the interplay between freshly 
injected relativistic particles and energy losses that cause a 
depletion of high-energy particles which results into a steep- 
ening of the spectrum at high frequencies. The discovery of 
a GPS radio source, considered to be a young object where 
the radio emission started a few thousand years ago, but dis- 
playing a steep spectral index (af;| > 1.0) in all its compo- 
nents is somewhat surprising. Indeed, the strong steepening 
found across the components of PKS 1518+047 is incompat- 
ible with the moderate high-frequency steepening predicted 
by a continuous injection of particles, suggesting that no in- 
jection/acceleration of fresh relativistic particles is currently 
occurring in any region of the source. Furthermore, when the 
injection of fresh particles is over, the strong adiabatic losses 
cause a fast shift of the spectral turnover towards lower fre- 
quencies and a decrement of the peak flux density. Therefore, 
detections of dying radio sources with a spectral peak still 
occurring in the GHz regime are extremely rare. A possible 
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Figure 3. VLA image at 1.6 GHz of PKS 1518+047. Another 
source, 1518+047B, is located at ~1 arcminute in the south-east 
direction from the target. 



explanation may arise from the presence of a dense ambient 
medium enshrouding the radio source which may confine 
the relativistic particles reducing the adiabatic losses. This 
is probably the case of fading radio sources found in galaxy 
clusters where the int racluster medium (ICM) is likely limit- 
ing adiabatic cooling IjSlee et al. 120011 ). However, in the case 
of PKS 1518+047 the source is surrounded by the interstel- 
lar medium (ISM) of the host galaxy, but its confinement 
cannot be related to neutral hydrogen, because observations 
searching for Hi absorption did not provide evidence of a 
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Table 2. Best fit parameters obtained with the CI and CI OFF models + SSA for the source PKS 1518+047 and for its northern and 
southern components. 



dense environment (IGupta et al. |[2006t ). 
From the analysis of the optically-thin spectrum of 
PKS 1518+047, we find that the steepening may be ex- 
plained assuming that the supply of relativistic plasma is 
still taking place, but the energy distribution of the injected 
particles is uncommonly steep. Another explanation is that 
the steep spectral shape is due to the absence of freshly in- 
jected/reaccelerated particles and the energy losses are driv- 
ing the spectrum evolution. Support to this scenario comes 
from the best fit to the spectrum of the northern compo- 
nent, where the continuous injection model fails in repro- 
ducing the spectral shape, even assuming that relativistic 
particles are injected with a steep energy distribution. This 
steep spectrum can be best reproduced by a synchrotron 
model in which no particle supply is taking place. Further- 
more, the time spent by the source in its "fader" phase is 
about 20% of the whole source age. 

We estimate the break energy 7b of the electron spectrum 
and the electron radiative time by: 

7 b~ 487^ /2 f^ 1/2 (l + Z ) 1/2 (2) 
and 

t s = 5.03 x W 4 H~ 3/2 u- 1/2 (l + Z y 1/2 (yr) (3) 
where H is in mG and v\, in GHz. 

If in Eqs. [2] and [3] we consider the equipartition magnetic 
field and v\, derived from the fits, we find that the energy 
of the electrons responsible for the break should be 7 ~ 
400 - 600 and the source radiative lifetime is t s — 2700 ± 
600 yr. As a consequence, electrons with 7 > 600 (i.e. with 
shorter radiative lifetime) have already depleted their energy 
and they contribute only to the high-frequency tail of the 
spectrum. 

If we consider the £off/£ s ~ 0.2, as derived from the fit 
to the overall spectrum, we find that £off = 550 ± 100 
years, and this represents the time elapsed since the last 
supply/acceleration of relativistic particles. These values 
indicate that the radio source was 2150±500 years old when 
the radio emission switched off. 

The small £off derived suggests that adiabatic losses had 
not enough time to shift the spectral peak of PKS 1518+047 
far from the GHz regime. In the presence of adiabatic 
expansion, and considering a magnetic field frozen in 
the plasma, the spectral peak is shifted towards lower 
frequencies: 



v P ,i = v pfi (—) (4) 

where u Pi o and v Pt i are the peak frequency at the time 
to and t\ respectively (|Orienti fc Dallacasa I l2008h . As 
the time elapsed after the switch off increases, the peak 
moves to lower and lower frequencies, making the source 
unrecognisable as a young GPS. For example, when toFF 
will represent more than 40% of the total source lifetime, 
the peak should be below 300 MHz, so it will be difficult to 
identify as a dying young radio source. Future instruments 
like LOFAR may unveil a population of such fading radio 
sources. 



5 CONCLUSIONS 

We presented results from multi-frequency VLBA and VLA 
observations of the GPS radio source PKS 1518+047. The 
analysis of the spectral index distribution across the whole 
source structure showed that all the source components 
have very steep optically-thin synchrotron spectra. The 
radio spectra are well explained by energy losses of rel- 
ativistic particles after the cessation of the injection of 
new plasma in the radio lobes. This result, together with 
the lack of the source core and active non-steep spectrum 
components like hot spots and knots in the jets, suggests 
that no injection/acceleration of fresh particle is currently 
occurring in any region of the source. For this reason, 
PKS 1518+047 can be considered a fading young radio 
source, in which the radio emission switched off shortly 
after its onset. When the supply of energy is switched off, 
given the high magnetic fields in the plasma, the spectral 
turnover moves rapidly towards low frequencies, making 
the source undetectable at the frequencies commonly used 
for radio surveys. However, in PKS 1518+047 the time 
elapsed since the last particle acceleration is of the order 
of a few hundred years, suggesting that PKS 1518+047 still 
has a GPS spectrum because adiabatic losses have not had 
enough time to affect the source spectrum, shifting the peak 
far from the GHz regime. If the interruption of the radio 
activity is a temporary phase and the radio emission from 
the central engine will re-start soon, it is possible that the 
source will appear again as a GPS without the severe steep- 
ening at high frequencies. If this does not happen, the fate 
of this radio source is to emit at lower and lower frequencies, 
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Figure 4. The best fit to the overall spectrum of PKS 1518+047 using the CI (left) and the CI OFF models (center), and the reduced 
chi-square versus the injection spectral index (right). Error bars are rather small, and often result within the symbol. 
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Figure 5. The best fit to the spectrum of the southern component Figure 6. The best fit to the spectrum of the northern component 

of PKS 1518+047 using the CI (top) and the CI OFF models of PKS 1518+047 using the CI (top) and the CI OFF models 

(bottom). (bottom). 
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until it disappears at frequencies well below the MHz regime. 
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